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Abstract 

The nonlinear Schrodinger equation possesses three distinct six-parameter families of complex- 
valued quasi-periodic travelling waves, one in the defocusing case and two in the focusing case. 
All these solutions have the property that their modulus is a periodic function oi x — ct for some 
c e M. In this paper we investigate the stability of the small amplitude travelling waves, both 
in the defocusing and the focusing case. Our first result shows that these waves are orbitally 
stable within the class of solutions which have the same period and the same Floquet exponent 
as the original wave. Next, we consider general bounded perturbations and focus on spectral 
stability. We show that the small amplitude travelling waves are stable in the defocusing case, 
but unstable in the focusing case. The instability is of side-band type, and therefore cannot be 
detected in the periodic set-up used for the analysis of orbital stability. 



Running head: Periodic waves in the NLS equation 

Corresponding author: Mariana Haragu§, haragus@math.univ-fcomte.fr 

Keywords: nonlinear Schrodinger equation, periodic waves, orbital stability, spectral stability 



1 Introduction 



We consider the one-dimensional cubic nonlinear Schrodinger equation (NLS) 

iUt{x,t) + U^^{x,t) ±\U{x,t)\'^U{x,t) = 0, 

where x G M, t £ M, C/(x, i) G C, and the signs + and — in the nonlinear term correspond 
to the focusing and the defocusing case, respectively. In both cases the NLS equation possesses 
quasi-periodic solutions of the general form 

U{x, t) = e'^P^-"^*) V{x-ct) , G M , t E M , (1.1) 

where p, w, c are real parameters and the wave profile F is a complex-valued periodic function of its 
argument. The aim of the present paper is to investigate the stability properties of these particular 
solutions, at least when the wave profile V is small. It turns out that the discussion is very similar 
in both cases, so for simplicity we restrict our presentation to the defocusing equation 

■iUt(x,t) + U^^{x,t) -\U{x,t)\^U{x,t) = 0, (1.2) 

and only discuss the differences which occur in the focusing case at the end of the paper. 

A crucial role in the stability analysis is played by the various symmetries of the NLS equation. 
The most important ones for our purposes are the four continuous symmetries: 

• phase invariance: U{x,t) i-^ U{x,t)e^'^, G M; 

• translation invariance: U{x,t) i— > U{x + ^,t), ^ G M; 

• Galilean invariance: + vt, t), V £R; 

• dilation invariance: U{x,t) XU{Xx, X'^t), A > 0; 

and the two discrete symmetries: 

• reflection symmetry: U{x,t) i-^ U{—x,t), 

• conjugation symmetry: U{x,t) U{x,—t). 

As is well-known, the Cauchy problem for equation (|1.2j) is globally well-posed on the whole 
real line in the Sobolev space Fi(M,C), see e.g. [HI lEl 1131 UHl • Alternatively, one can solve the 
NLS equation on a bounded interval [0, L] with periodic boundary conditions, in which case an 
appropriate function space is ifpgj,([0, L], C). In both situations, we have the following conserved 
quantities: 

EiiU) = ^l^\U{x,t)\^dx , 

E2{U) = ^J^U{x,t)U^{x,t)dx, 

EsiU) = J^(^^\U,ix,t)\^ + \\Uix,tt)dx , 
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where X denotes either the whole real line or the bounded interval [0, L]. The quantities Ei and 
E2 are conserved due to the phase invariance and the translation invariance, respectively, whereas 
the conservation of £^3 originates in the fact that equation ()1.2() is autonomous. 

The symmetries listed above are also useful to understand the structure of the set of all quasi- 
periodic solutions of ()1.2I) . Assume that U{x,t) is a solution of (|1.2() of the form where 
1/ : M — > C is a bounded function. Since \U{x,t)\ = \V{x — ct)\, the translation speed c € M is 
uniquely determined by U, except if the modulus \ V\ is constant. In any case, using the Galilean 
invariance, we can transform U {x, t) into another solution of the form 1)1.1(1 with c = 0. Once this 
is done, the temporal frequency lo is in turn uniquely determined by U{x,t), except in the trivial 
case where V is identically zero. In view of the dilation invariance, only the sign of u is important, 
so we can assume without loss of generality that lo £ {— 1;0; 1}. Setting U{x,t) = e~"^*iy(x), we 
see that W{x) = e^^^V{x) is a bounded solution of the ordinary differential equation 

W^^{x) + LoW{x) - \W{x)\'^W{x) = , xeR. (1.3) 

If a; = or (J = —1, is it straightforward to verify that W = is the only bounded solution of 
TT^ . thus we assume from now on that lo = 1. Equation (|1.3|) is actually the stationary Ginzburg- 
Landau equation and the set of its bounded solutions is well-known IHl El • There are two 
kinds of solutions of (jl.3j) which lead to quasi-periodic solutions of the NLS equation of the form 
CU): 

• A family of periodic solutions with constant modulus W{x) = (1 — p^'j^/'^ e^(P^+'^'\ where 
p G [—1, 1] and ip £ [0, In]. The corresponding solutions of (|1.2() are called plane waves. The 
general form of these waves is 

u{x, t) = S^^-^^) V , 

where p G M, G M, and V £ C satisfy the dispersion relation 00 = p^ -\- \ V\^ . 

• A family of quasi-periodic solutions of the form W{x) = r(x)e"^^^-', where the modulus r(x) 
and the derivative of the phase ip{x) are periodic with the same period. Any such solution can 
be written in the equivalent form W{x) = e'^^ Q{2kx), where p G M, A; > 0, and Q : M ^ C 
is 27r-periodic. In particular, 

U{x, t) = e-''W{x) = e'(P^-*)Q(2A:x) (1.4) 

is a quasi-periodic solution of (|1.2j) of the form ()1.1|) (with c = and lo = 1). We shall refer 
to such a solution as a periodic wave, because its profile \U{x,t)\ is a (non-trivial) periodic 
function of the space variable x. Important quantities related to the periodic wave (jl.4j) are 
the period of the modulus T = Tr/k, and the Floquet multiplier e^^"^. For small amplitude 
solutions {\Q\ <^ 1) the minimal period T is close to vr, hence A; ~ 1, and the Floquet 
multiplier is close to —1, so that we can choose p ~ 1. 

While the plane waves form a three-parameter family, we will see in Section |21 that the periodic 
waves form a six-parameter family of solutions of (|1.2() . However, the number of independent 
parameters can be substantially reduced if we use the four continuous symmetries listed above. 
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Indeed it is easy to verify that any plane wave is equivalent either to Ui{x,t) = or to U2{x,t) = 
e~'*. In a similar way, the set of all periodic waves reduces to a two-parameter family. 

As far as the stability of the plane waves is concerned, the conserved quantities Ei and E3 
immediately show that the trivial solution Ui = is stable (in the sense of Lyapunov) with respect 
to perturbations in H^(M.) or i?pgj.([0, L]), for any L > 0. The same conservation laws also imply 
that the plane wave U2 = e~'* is orhitally stable in the following sense. Assume that I = [0, L] is 
a bounded interval, and let U{x,t) be the solution of ()1.2|) with initial data U{x,0) = 1 + Vo(x), 
where Vq G H^^^iX) and || Vo||j:fi(x) < e. If e > is small enough, then 

inf \\U{-,t)-e"^\\Him < C{l)e , for ah t G M , (1.5) 

V5e[0,27r] 

where the constant C{I) depends only on the length of the interval X. This stability property is 
easily established using the conserved quantity 

E{U) = jj^\\U,{x,t)\'' + \{\U{x,tt -lf)dx = E:,{U) - E^{U) + . 

A similar result holds for small perturbations of U2 in H^{W). In that case, the bound H1.5|) 
holds for any bounded interval X C M, but the conservation of E(U) does not prevent the norm 
||[/(-,t) — e~'*||j|/i(iR) from growing as \t\ 00. We refer to |2H1 Section 3.3] for a detailed analysis 
of the stability of plane waves. 

The stability question is much more difficult for periodic waves. In contrast to dissipative 
systems for which nonlinear stability of periodic patterns has been established for rather general 
classes of perturbations, including localized ones (see e.g. HHD) no such result is available so far for 
dispersive equations. In the particular case of NLS, the stability of the ground state solitary waves 
has been intensively studied (see e.g. [3123)) but relatively little seems to be known about the 
corresponding question for periodic waves. A partial spectral analysis is carried out by Rowlands 
|23j , who shows that the periodic waves are unstable in the focusing case and stable in the defocusing 
case, provided disturbances lie in the long-wave regime. Spectral stability has also been addressed 
for certain NLS-type equations with periodic potentials [Hll^. Very recently, Angulo P has shown 
that the family of "dnoidal waves" of the focusing NLS equation is orbitally stable with respect to 
perturbations which have the same period as the wave itself. In all these previous works, the wave 
profile V is assumed to be real-valued. Here we restrict ourselves to small amplitude solutions, 
but allow for general complex- valued wave profiles. While the nonlinear stability of these waves 
with respect to bounded or localized perturbations remains an open problem, we treat here two 
particular questions: orbital stability with respect to periodic perturbations, and spectral stability 
with respect to bounded or localized perturbations. 

Our first result shows that the periodic waves of H1.2|) are orbitally stable within the class of 
solutions which have the same period and the same Floquet multiplier as the original wave: 

Theorem 1 (Orbital stability) 

Let X = i/ppj.([0, 27r], C). There exist positive constants Co, eo, and 60 such that the following 
holds. Assume that W{x) = e'^^ Qpcr(2A;x) is a solution of Ill.^A) with Qpcr G X, ||Qper||x < Sq, 
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and p,k « 1. For all R ^ X such that \\R\\x < co? ihe solution U{x,t) = e'*^^^ *^Q{2kx,t) of M.'J(l 
with initial data U{x, 0) = e'^^(Qper(2/cx) + R{2kx)) satisfies, for all t G M, 

inf ||Q(.,t)-e''^Qper(--e)IU < Copllx . (1-6) 

■,5,?e[0,27r] 



Remarks 

1. We point out that Theorem ^ holds uniformly for all quasi-periodic solutions of (|1.2() with small 
amplitude. In particular the unperturbed solution e'^^'^"*) Qperi^kx) can be either a periodic wave, 
or a plane wave, or even the zero solution. 

2. The proof of Theorem ^ relies on the classical approach to orbital stability which goes back to 
Benjamin ^ (see also |2 El El)- While for solitary waves this method gives a rather complete 
answer to the stability question, in the case of periodic waves it allows to prove orbital stability 
only if we restrict ourselves to solutions which have the same periodicity properties as the original 
wave (see however Remark 13 . 1 1 1 b elow for a discussion of this limitation). In this paper we shall 
use the general framework developed by Grillakis, Shatah, and Strauss ^lEl) with appropriate 
modifications to obtain a uniform stability result for small waves. Note that a direct application 
of the stability theorem in would give the same conclusion as in Theorem ^ but with stability 
constants Cq and eo depending on the wave profile Qper- 

3. Following the approach of ^21 it is shown in jT^ that all periodic waves of ()1.2() are orbitally 
stable in the sense of (|1.6|) . without any restriction on the amplitude of the wave profile Qper- 
The argument in ^] relies in part on the results obtained in the present paper, and uses a global 
parametrization of the set of quasi-periodic solutions of (|1.3j) which is very different from the 
explicit series expansions that we use here to describe the small amplitude solutions. 

4. It is worth considering what Theorem ^ exactly means in the particular case where TV is a 
real-valued periodic solution of (|1.3j) (such a solution is often referred to as a "cnoidal wave" in 
the literature). In that case we have W{x) = e^^^ Qperi^kx) where p = k = n/T and T > vr is 
the minimal period of \W\. The Floquet multiplier e^^"^ is therefore exactly equal to —1, so that 
W{x + T) = —W{x) for all x G M. In particular, we see that W is periodic with (minimal) period 
L = 2T. Thus Theorem n shows that the L-periodic cnoidal wave U{x,t) = e~'*iy(x) is orbitally 
stable with respect to L-periodic perturbations W provided that W{x + L/2) = —W{x) for all 
X G R. As is explained in [T|, without this additional assumption the classical approach does not 
allow to prove the stability of cnoidal waves with respect to periodic perturbations. 

Next, we investigate the spectral stability of the periodic waves with respect to bounded, or 
localized, perturbations. Although spectral stability is weaker than nonlinear stability, it provides 
valuable information about the the linearization of the system at the periodic wave. Our second 
result is: 

Theorem 2 (Spectral stability) 

Let Y = L^(IR, C) orY = Cfe(M, C). There exists 5i > such that the following holds. Assume that 
W{x) = e^P^ Qper{2kx) is a solution of il.,'^) with Qper £ X, ||<5por||x < Si, and p,k ^ 1, just as 
in Theorem^ Then the spectrum of the linearization of hl.^) about the periodic wave e~'*VF(x) in 
the space Y entirely lies on the imaginary axis. Consequently, this wave is spectrally stable in Y . 
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The proof of Theorem |2l is based on the so-called Bloch-wave decomposition, which reduces the 
spectral study of the linearized operator in the space Y to the study of the spectra of a family 
of linear operators in a space of periodic functions. Bloch waves are well-known for Schrodinger 
operators with periodic potentials H2] and have been extensively used in dissipative problems 
pni OH 051 OB] , but also in a number of dispersive problems [HI El 1^- The advantage of such 
a decomposition is that the resulting operators have compact resolvent, and therefore only point 
spectra. The main step in the analysis consists in locating these point spectra. For our problem, we 
rely on perturbation arguments for linear operators in which we regard the operators resulting from 
the Bloch-wave decomposition as small perturbations of operators with constant coefficients. The 
latter ones are actually obtained from the linearization of (|1.2|) about zero, and Fourier analysis 
allows to compute their spectra explicitly. The restriction to small amplitudes is essential in this 
perturbation argument, and we do not know whether spectral stability holds for large waves. 

The rest of the paper is organized as follows. In Section |2l we briefly describe the set of all 
bounded solutions of (|1.3|) . and we introduce an analytic parametrization of the small amplitude 
solutions which will be used throughout the paper. In Section El we recall the main ideas of the 
orbital stability method, and we apply it with appropriate modifications to prove Theorem ^ 
Spectral stability is established in Section |11 using Bloch-wave decomposition and the pertubation 
argument mentioned above. 

Finally, in Section |S1 we discuss the stability of the small periodic waves of the focusing NLS 
equation. In contrast to the defocusing case, the focusing NLS equation possesses two different 
families of quasi-periodic solutions of the form one for u; > and the other for oj < 

Small solutions exist only within the first family, and their stability properties can be analyzed 
as in the defocusing case. However, while for periodic perturbations we obtain the same orbital 
stability result as in Theorem Q it turns out that the small periodic waves are spectrally unstable in 
the focusing case. Unstable spectrum is detected for perturbations with wave-numbers which are 
close to that of the original wave (side-band instability) . As for the second family, which contains 
only large waves, we refer to ^ for a proof of orbital stability and to [2] for a discussion of 
spectral stability. 

Acknowledgments The authors thank A. De Bouard and L. Di Menza for fruitful discussions. 
This work was partially supported by the French Ministry of Research through grant ACI JC 1039. 

2 Parametrization of small periodic waves 

In this section, we briefly review the bounded solutions of equation H1.3|) with lo = 1: 

W^^{x) +W{x) -\W{x)\'^W{x) = , W:R^C, (2.1) 

and we give a convenient parametrization of all small solutions. If we interpret the spatial variable 
X G M as a "time" , equation (|2.1|) becomes an integrable Hamiltonian dynamical system with two 
degrees of freedom. The conserved quantities are the "angular momentum" J and the "energy" E: 

J = Im(W,) , E = i|iy,|2 + _ ^\W\^ . (2.2) 
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If is a solution of t\2.1\i with J ^ 0, then W{x) ^ for all x G M, so that we can introduce the 
polar coordinates W{x) = r{x) e^'^^^\ The invariants then become 

9 ^ ,P 
J = r-'m E = —-\ . 



The set of bounded solutions of (|2.1|) can be entirely described in terms of these two invariants 
[HlElIinilll]- In the parameter space (J, E) there is an open set 

D = {(J,^) G j2 < 4/27, E^{J) <E < S+(J)} , (2.3) 

where E-,E^ are explicit functions of J, such that the closure D consists of all values of {J,E) 
which give rise to bounded solutions W of (|2.H) (see Fig. 1). Furthermore, we have the following 
classification for {J,E) in D: 

(i) If = E^{J), then is a periodic solution with constant modulus and linear phase, i.e. 
W{x) = Wp^^ix) = (1 _p2)i/2ei{p^'+'^) with 1/3 < < 1 and ip £ [0,2tt]. 

(ii) If E = Ej^{J), then either W = Wp^^ for some < 1/3 and some ip E [0, 27r], or W is a 
homoclinic orbit connecting Wp^^_ at x = — oo to Wp^^^ at x = +oo for some ip-., G [0, 27r]. 

(iii) If E-{J) < E < E^{J) and J ^ 0, then the modulus and the phase derivative of W are 
both periodic with the same period T(J,E) > vr. Let ^{J,E) be the increment of the phase 
over a period of the modulus, so that W{x + T) = e^^W{x) for all x G M. In general $ is 
not a rational multiple of tt, hence the solution W of (|2.1j) is typically not periodic, but only 
quasi-periodic. In the particular case where J = 0, then e'* = — 1 and W is periodic with 
period 2T{0,E). 



1 
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Fig. 1: The region D in the parameter space {J,E) for which H2.1|l has bounded solutions. 

For a fixed pair (J, E) G D, the bounded solution W of (|2.1|) satisfying (|2.2jl is unique up to a 
translation and a phase factor. In case (iii), the period T and the phase increment $ (or the Floquet 
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multiplier e'*) are important quantities which play a crucial role in the stability analysis of the 
quasi-periodic solutions of ()2.1() . both for the Schrodinger and the Ginzburg-Landau dynamics. A 
number of properties of T and $ are collected in T!T^. In particular, if we define the renormalized 
phase 

I if J = , 

then T : D ^ M and ^ : D ^ M are smooth functions {J, E) £ D, in contrast to ^( J, E) which 
is discontinuous at J = 0. In addition, T ~ vr and \I' ~ for small solutions ~ 0. 



The periodic solutions Wp^^, of (|2.1|) correspond to plane waves of the NLS equation. We are 
mainly interested here in the quasi-periodic solutions described in (iii) above, which correspond to 
periodic waves of the form (|1.H) . To see this, fix {J,E) G D and let PV^ : M — > C be a bounded 
solution of (|TT|) satisfying We set 

W{x) = e'^^P(A:x) , xeR, (2.5) 

in which k and £ are related to the period T{J,E) and the renormalized phase ^{J,E) through 

As |Ty(x)| = \P{kx)\ is T-periodic (in x) by the definition of T{J,E), it is clear that \P{y)\ is tt- 
periodic (in y). Moreover, since W{x+T) = e^^W{x) = —e^'^W{x), we also have P(y + 7r) = —P{y) 
for all y G M, hence P is 27r-periodic. Thus U{x, t) = e~^^W{x) = e^^^e~^^P{kx) is a quasi-periodic 
solution of (|1.2j) of the form with = 1 and c = 0. 

Remark 2.1 Using the continuous symmetries of the NLS equation, we obtain for each pair 
(J, E) G D a four-parameter family of periodic waves: 

Uc,x,^,l^{x,t) = Ae'P->^"e-'"-^*e'^P(A;A(x-ct)-e) , 

where c € M, A > 0, (/?, ^ G [0, 27r], and Pc^\ = Xi + c/2, w^a = )? + cXl + c^ / A. Taking into account 
the parameters J,E, we obtain altogether a six-parameter family of periodic waves of hl.^) . 



Alternatively, we can write the solution (|2.5)) of (|2.1j) in the form 

W{x) = e^(^+^)^Q+(2A:x) = e'^^''^)^ g-(2A;a;) , x G M , (2.7) 

where Q^{z) = e^'^/^P(z/2). By construction, and |Q^| are now periodic functions with the 
same minimal period 27r. The representation (|2.7|) turns out to be more convenient than (|2.5)) to 
study the orbital stability of the periodic waves in the next section. 

The global parametrization of the quasiperiodic solutions of (|2.1|) in terms of the invariants 
(J, E) is natural, but it is not very convenient as far as small solutions are concerned because the 
admissible domain D is not smooth near the origin (see Fig. 1). For this reason, we now introduce 
an analytic parametrization of the small solutions of (|2.1j) . We start from the representation (|2.5|1 , 
and we choose as parameters the first nonzero Fourier coefficients of the 27r-periodic function P: 

"^^C ^^^^ ' ^ " i C ^^^^ ""^'^ ■ 
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(Remark that P has zero mean over a period.) Replacing P{y) with e^"^P(y + .^) if needed, we can 
assume that both a and b are real. If P (hence also W) is small, we have T ~ vr and ^' ~ 0, hence 
/c ~ 1 and ^ ~ 0. This determines uniquely the expansion of P, fc, H. in powers of a and h. Setting 

Wa,b{x) = e'''^'^^Pa,bika,bX) , xGR, (2.8) 

we obtain after straightforward calculations: 

ia,b = l{b^-a^) + 0{a^ + b^) , 

ka,b = l-^(a2 + 62) + o(a^ + 6^) , (2.9) 

PaM = ae-'y + be'y - —e-^'y -—e'^'y + 0{\ab\{\a\' + \b\')) , 

o o 

as {a,b) — > (0,0). Notice also that the invariants J,E have the following expressions: 

J = b''-a^ + ]^{a^-b^) + 0{a^ + b^) , 

E = a2 + 62_3aV-^(a^ + &^) + 0(a6 + 66) . 

With this parametrization, replacing a with —a or 6 with —b gives the same function P up to a 
translation and a phase factor: 

P-aM = -'^Pa,b{y + Vr/2) , P.a,-b{y) = 'PaM = Pa,b{y + Vt) , y G M . 

It follows that J, -E, hence also A;,£, are even functions of a and b. Similarly, Pb,aiy) = Pa,b{y)- 
This conjugation leaves E unchanged but reverses the sign of J, hence ka^ = kb^a and la,b = —h,a- 
Therefore, using the symmetries of 1)2. we can restrict ourselves to the parameter region {b > 
o > 0} without loss of generality. 

Two particular cases will play a special role in what follows. 

(i) (Cnoidal waves) If a = 6, then ia,a = and we obtain a family of real-valued periodic 
solutions Wa,a{x) = Pa,a{ka,ax), where 

ka,a = l-^a2 + 0(a^), Pa,a{y) = 2a COS y - ^ cos(3y) + 0(|a|^) . 
Observe that J = in that case. 

(ii) (Plane waves) If a = 0, then Pofi{y) = be^y, hence WQ^b has constant modulus. It follows 
that 

Wo,b{x) = be'-^^ , and ko,b + lo,b = ^Jl-b^ ■ 
In addition, one also finds k^^b = \/l — 36^/2. It is advantageous here to use the representa- 
tion dSIl), namely W^^biy) = e'P^''^Q^fe(y) with = £o,6 + ko,b = (1 - &')'/' and = b. 
Similarly, if 5 = we have Pa,oiy) = oe~'^ and thus 

Wa,oix) = ae-'^'^" , and ka,o - 4,o = Vl-a^ , = \/l - 3aV2 ■ 

Alternatively, Wafi{y) = e^'^'^yQafiiv) with p'^ = 4,o - ka,o = -(1 - a^^^^ and Q^q = a. 
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3 Orbital stability 



In this section we prove the orbital stabiUty result in Theorem ^ Since we restrict ourselves to 
periodic waves with small amplitude, we shall use the local parametrization (|2.8() . ()2.9|) of the small 
solutions of ()2.1|) . Given (a, 6) G with ||(a,6)|| sufficiently small, we consider the periodic wave 
Ua,b{x,t) = e-'^Wa,b{x), where 

Wa,b{x) = e''^'>'^Pa,b{ka,bX) = e'^^'"^ Q a,b{2ka,bx) , X G R . 

Here ia,b,ka,b,Pa,b are defined in ()2.9() . and the last expression in the right-hand side corresponds 
to the first choice in (|2.7p . namely 

Pa,b = ia,b + ka,b , Qa,b{^) = (T'^'^ Pa,b{z /2) . (3.1) 

From the properties of Pa b we deduce 



Q^a,b{z) = Qa,b{z + Tt), Q_a-b{z) = -Qa,b{z), Qb,a{^) = e '''Qa,b{z) , (3.2) 

and that the real and imaginary parts of Qa,b are even and odd functions of z, respectively. 

Remark 3.1 Without loss of generality, we shall assume henceforth that 6 > a > 0. Note that the 
second choice in \2. ?| ) would he preferable when >lP'. 

3.1 Main result and strategy of proof 

To study the stability of Ua,b{x,t) we consider solutions of 1)1.2(1 of the form 

U{x,t) = e^(P-''^-*)Q(2fca,ba;,t) , (3.3) 
where Q{z,t) is a 27r-periodic function of z which satisfies the evolution equation 

iQt + ^iPa,bka,bQz + ^klbQzz + (1 - plb)Q - \Q?Q = • (3.4) 

By construction, Qa,b{z) is now a stationary solution of (|3.4|) and our goal is to show that this 
equilibrium is stable with respect to 27r-periodic perturbations. We thus introduce the function 
space 

X = Fji<,,([0,27r],C) = |n G H^^^i^X) n(z+27r) = u{z) , Vz G m} , 
which is viewed as a real Hilbert space equipped with the scalar product 

{u,v)x = Re / {u{z)v{z) + Uz{z)vz{z))dz , u,v ^ X . 







As usual, the dual space X* will be identified with F~i,([0, 27r], C) through the pairmg 



27r 



(n, v) = Re u{z)v{z) dz , u G X* , v X . 
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It is well-known that the Cauchy problem for H3.4|) is globally well-posed in the space X. More- 
over, the evolution defined by (|3.4|) on X is invariant under a two-parameter group of isometries. 
The symmetry group is the two-dimensional torus G = = (IR/27rZ)^, a compact abelian Lie 
group which acts on X through the unitary representation IZ defined by 

(7^(^,^)n)(z) = e--^^u{z + i) , ueX , {ip,O^G. 

Due to these symmetries, it is useful to introduce the semi-distance p on X defined by 

p{u,v) = inf \\u — TZr^ ^)v\\x , u,v G X . (3-5) 

(</p,C)eG ^^'^^ 

In words, p(n, v) is small if u is close to v (in the topology of X) up to a translation and a phase 
rotation. Our stability result in Theorem^ can now be formulated as follows: 

Proposition 3.2 There exist Cq > 0, eo > 0, and 6o > such that, for all (a, b) G M? with 
||(a,6)|| < 6o, the following holds. If Qq G X satisfies p{Qo,Qa,b) < e for some e < eo, then the 
solution Q{z, t) of \'3.4^ with initial data Qq satisfies p{Q{-,t), Qa,b) ^ Cqe for all t G R. 

For each fixed value of (a, 6), the stability of the periodic (or plane) wave Qa^b can be proved 
using the abstract results of Grillakis, Shatah and Strauss jl4l I15j . However, this approach would 
not give a stability theorem that holds uniformly in a neighborhood of the origin, as it is the case in 
ProDosition l3.2l A difficulty in proving such a uniform result is that we have to deal simultaneously 
with three sorts of solutions: the zero solution (a = 6 = 0), plane waves {ah = 0) and periodic 
waves [ah ^ 0). These equilibria are genuinely different from the point of view of orbital stability 
theory, because their orbits under the action of the symmetry group G have different dimensions 
(0, 1, and 2, respectively). In what follows, we shall concentrate on the periodic waves, and at the 
end we shall indicate how the other cases can be incorporated to obtain a uniform result. Whenever 
possible, we shall adopt similar notations as in to facilitate comparison. 

Due to its symmetries, equation 1)3. 4|) has the same conserved quantities as the original NLS 
equation, namely 

M{Q) = rQ{z)Q,{z)dz , 
^ Jo 

£{Q) = l^'\2kl,mz)\' + \\Qizt)dz. 

The charge A^, the momentum M and the energy £ are smooth real-valued functions on X. Their 
first order derivatives are therefore smooth maps from X into X*: 

N'{Q) = Q , M'iQ) = iQ, , £'{Q) = -Akl^Q,, + IQI^Q . 

Similarly, the second order derivatives are smooth maps from X into C{X,X*), the space of all 
bounded linear operators from X into X*: 

N"{Q) = 1 , M"(Q) = id, , £"{Q) = -Akl.d,, + \Q\^ + 2Q ® Q , 
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where 

f27r 



{{Q®Q)u,v) 



/ Re(Qli) Re{Qv) dz, y u,v e X. 
Jo 



From now on, we fix {a,b) £ M? with ||(a,6)|| sufficiently small. As is explained above, we 
assume for the moment that a6 7^ 0, in which case the function Qa,b £ ^ defined by is a 

stationary solution of H3.4() corresponding to a periodic wave of the original NLS equation, i.e. 
\Qa,b\ is not constant. By construction, Qa^b is a critical point of the modified energy 

£a,b{Q) = £iQ) - (1 - Plb)N{Q) - Apa,bka,bM{Q) , (3.6) 

namely <?^^((5a,b) = 0. The orbital stability argument is based on two essential ingredients: 

Claim 1: The equilibrium Qa^b is a local minimum of the function 8a,h restricted to the codimension 
two submanifold 

Sa,f>(Q) = {q G X I N{Q) = N{Qa,b) , M{Q) = M{Qa,b)] • (3.7) 

Note that this manifold contains the entire orbit of Qa,b under the action of G. 

Claim 2: The equilibrium Qa,b is a member of a two-parameter family of travelling and rotating 
waves of the form 

Q(z, t) = e"^^* Q'^^ll^z + ct) , zGR, t£R, (3.8) 

where (w, c) lies in a neighborhood of the origin in M? (the Lie algebra of G) and Q^'^ G X is a 
smooth function of {u!,c) with = Qa,b- Moreover the map {u!,c) {Q'^'^) , M {Q^'^)) is a 

local diffeomorphism near (w, c) = (0,0). 



3.2 Proof of Claim 2 

The second claim is easily justified using the continuous symmetries of the NLS equation. Indeed, 
let {a',b') G be close to (a, 6). Then 

U(x,t) = SPa'M^-^) Qa',b'{^ka',b'x) , X G M , t G M , 

is a solution of the NLS equation, but it is not of the form 1)3. 3|) because pa',b' 7^ Pa,b and ka'^b' 7^ kafi 
in general. However we can transform U{x,t) into a solution of H1.2|) of the form 1)3. 3() . 1)3. 8|) by 
applying successively a dilation of factor A and a Galilean transformation of speed v, where 

^ \a.',b' kab a'.b' r,/\a',b' \ n\ 

^ = K b = 1 — , V = v^ ^^ = 2(A 'j, pa',b' - Pa,b) ■ (3.9) 

'^a',b' 

After some elementary algebra, we obtain Q^'l{z) = A^'^ Qa',b'{z) with 



^ = - P<V) - (1 - Pa,fc) , C = 4(A;',' )'A:a',fe'Pa',6' " ^ka,bPa,b ■ (3-10) 

Using the expansions (|2.9() . it is straightforward to verify that 



duj 


dc 


IW 


'UaJ 


duj 


dc 


W 


W 



{a',b')={a,b) 
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4a —2a\ , o ,o 



26 26 .>1+C(a^+'''))- 



Since we assumed that ab ^ 0, the matrix Aia,b is invertible for a, b sufficiently small, hence the 
mapping {a',b') i-^ (w,c) defined by (|3.1()|) is a diffeomorphism from a neighborhood of (a, 6) onto 
a neighborhood of (0,0). This proves the existence of the travelling and rotating wave (|3.8|) for 
(u), c) S M? sufficiently small. Remark that the proffie Q^'^ is a critical point of the functional 

£a:biQ) = £a,biQ) - ojN{Q) - cM{Q) , QeX . 

Following ^H]) we define da,b{^,c) = ^a biQab)- "^^^ properties of the function da^b will play an 
important role in the orbital stability argument. 



Lemma 3.3 The Hessian matrix of the function da^b satisfies: 

n 



def / diu''' duidc 



dcduj dc'' 



(a;,c)=(0,0) 



TT / -2 -1 

3 I -1 1 



(l + 0(a2+62)) 



Proof. Since is a critical point of iS^ ^ , we have 



^da,b{u^,c) = -iV(Q^p , §-Ja,b{^,c) = -M{Qll) . (3.11) 

To compute the second-order derivatives, we parametrize {uj,c) by {a',b') as above. Using 1)3. 10|) 
we find 7ia,b = -iMa,by^lCa,b, where 

(irmZb) i^MiQ^fb) 



a,b 



{a' ,b')={a,b) 



Qab = Kb Qa',b'{z), we have 



mZb) = i^ibfma',b') , m{q:;^) = (a2Vm(q,,,o . 

On the other hand, using the expansion 

Qa,b{z) = ae-^^ + 6-4^e-2--^e- + 0(|a6|(|a|3 + |6|3)) , (3.12) 

o o 

which follows from ()2.9() . 1)3. 1(1 . we easily find 

NiQa,b) = vr(a2 + 62) + 0{aH\a^ + 6^)) , M(Q^,b) = vra^ + 0{aH\a^ + b^)) . 
Combining these results, we obtain 

(^■■'''"' = 6^(-\ ;)(i + o("^+''^)). >c.,-2.[i ;)^^^o^a^W)), 

and the conclusion follows. □ 

Lemma 13.31 implies that the Hessian matrix 'Ha,b is nondegenerate for ||(a,6)|| sufficiently small 
(in fact, 7ia,b has one positive and one negative eigenvalue.) It follows that the map {u:,c) ^ 
(A^((5^'^), M((5^'^)) is a local diffeomorphism near {uj,c) = (0,0), because by (|3.11|) the Jacobian 
matrix of this map at the origin is just —Ti.a,b- Thus Claim 2 above is completely justified. 
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Remark 3.4 At this point we could apply the general result of 115]/ . but as already mentioned this 
would not give the uniform result in Theorem According to the Stability Theorem in 11 Sf . in 
order to establish the orbital stability of a single wave Qa,b with ab ^ it suffices to show that the 
linear operator 

Ha,b = £abiQa,b) = -^k^ ^8,, - Aipa,bka,bdz - {1-pl J + \Qa,b\^ + 2Qa,b ^ Qa,b , (3.13) 



u d 

^^(V,5)<3a,fe , = -lQa,6 , -^T^{ip,i)Qa,b ^^^^^ = ^^^^^^ , (3.14) 



has precisely one simple negative eigenvalue, a two-dimensional kernel spanned by 
d_ 

and that the rest of its spectrum is strictly positive. Observe that Ha^b is self-adjoint in the real 
Hilbert space Lppj.([0, 27r], C) equipped with the scalar product (•,•). Clearly, the vectors \3.14\! 
always belong to the kernel of Ha^b- In fact, for small {a,b) we can determine the spectrum of Hafi 
by a perturbation argument similar to the one used for the spectral analysis of the operators Aa,b,'y 
in Section ^ We find that Hafi has exactly four eigenvalues in a neighborhood of the origin, the 
rest of the spectrum being positive and bounded away from zero. Among these four eigenvalues, 
two are always zero, and the other two have negative product —12a^b'^{l-\-0{a^-\-b'^)). This implies 
that Ha^b has the required properties, so that the wave profile Qa,b is orbitally stable if ab ^ 0. This 
information on the spectrum of Ha^b will not be used in the remainder of this section. However, 
since it provides the starting point for the stability analysis of large waves in ^11], we give a brief 
proof in the Appendix. 

3.3 Proof of Claim 1 

We now turn back to Claim 1 and study the behavior of the energy £a,b on the manifold T,afi 
defined by p.7|) . In the arguments below, we assume 6 > a > 0, so exclude for the moment the 
plane wave corresponding to a = 0. Let Ta^b be the tangent space to Sa^f, at the point Qa,b- 

ra,b = {Qex\ {N'{Qa,b),Q) = {M'{Qa,b),Q) = o} . 

Then X = Ta^b © J^a,b-, where J\fa,b (the "normal" space) is the two-dimensional subspace of X 
spanned by N'{Qa,b) = Qa,b and M'{Qa,b) = idzQa,b- When (a, b) is small, a more convenient basis 
of J\fa,b is {L,b,Va,b}, where 

= -dzQa,b = e-'' + 0{\ab\ + 62) , rja,b = ^(Qa,6 - iS.Qa.b) = 1 + 0{a^ + |a6|) . (3.15) 
a b 

The tangent space is further decomposed as Ta^b = Ya,b © Za^b, where 

Ya,b = {Q G Ta,b I {lQa,b, Q) = {dzQa,b, Q) = o} , 

and Za^b is the two-dimensional space spanned by iQa,b and dzQa,b- In view of (|3.14|) . Za^b is just 
the tangent space to the orbit of Qa,b under the action of G. Again, a convenient basis of Za^b is 
{i^a,b,i??a,b}. 

As in ^15], we introduce an appropriate coordinate system in a neighborhood of the orbit of 
Qa b under the action of G: 
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Lemma 3.5 Assume that ||(a,6)|| is sufficiently small and b > a > 0. There exist k > 0, Ci > 0, 
and C2 > such that any Q £ X with p{Q, Qa,b) < 1^0. can be represented as 

Q = T^(^,oiQa,b + j^ + y) , (3.16) 

where (ip,^) £ G, u £ Ma,b, V S and \\iy\\x + \\y\\x < Cip{Q,Qa,b)- Moreover, if Q £ T.a,b, 
then \\iy\\x < {C2/a)\\y\\j^. 

Remark 3.6 Here and in the sequel, all constants Ci,C2,... are independent of (a, 6) provided 
II (a, 6) II is sufficiently small. 

Proof. It is clearly sufficient to prove the result for all Q G X with ||Q — Qa, fe||x ^ ^a, where 
K > is a (small) constant that will be fixed below. Since X = Ma,b © ^a,fe © Za,b, any such Q can 
be decomposed as 

Q = Qa,b - XliQa,b + X2dzQa,b + 1^1 + Vl , 

where i^i £ Ma^b-, Hi G ^a,6; and (xi,X2) £ is the solution of the linear system 

(iQa,6, iQa,6) — ^^Qa,b, dzQa,b) \ [ Xl\ _ f — {^Qa,b, Q — Qa,l 



(3.17) 

-{dzQa,b,^Qa,b) {dzQa,b, dzQa,b) J J \ {dzQa,b, Q — Qa,b) J 

The matrix Aia,b in the left-hand side of (|3.17|) is invertible, and using the expansions (|.S.12|) we 

Since 6 > a, it follows that |j;i| + |x2| < {C/a)\\Q — Qa,b\\x < Ck for some C > (independent of 
o, 6). Now, since 

T^(ip,OQa,b = Qa,b - ^'^Qa,b + ^dzQa,b + 0{ip^ + S,"^) , 

the Implicit Function Theorem implies that, if {xi,X2) is sufficiently small, there exists a unique 
pair {ip,£,) £ ^ith = (xi,X2) + 0{xl + xj) such that T^^^j:)Q - Qa,b £ J^a,b © ^a,fe (see 

Lemma 4.2 in ^S] for a similar argument). Setting TZ'^^^-^Q — Qafi = i' + y, we obtain the desired 
decomposition (assuming that k > is small enough so that we can apply the Implicit Function 
Theorem). This choice of (92,^) does not minimize the distance in X between Q and TZ(^^^^-^Qa,b, 
because the subspaces Ma,b, Zafi, and Ya^b are not mutually orthogonal for the scalar product of 
X. However, since the minimum gap between these spaces is strictly positive (uniformly in a,b), 
we still have ||z^||x + ||y||x < C'll'^ + y\\x < Cip{Q,Qa,b)- (We refer to JS| for the definition and 
the properties of the minimum gap between closed subspaces of a Banach space.) 

Now, we assume in addition that Q £ Sa,fe) i.e. N{Q) = N{Qafi) and M{Q) = M{Qa,b)- In 
view of (|3.1(i)) . we have 

N{Q) = N{Qa,b + i^ + y) = N{Qa,b) + {Qa,b, V + y) + N {u + y) , 

and using the fact that y £ Ya^b C Ta^b we obtain {Qa,b,i^) + ^{^^ + y) = 0. A similar argument 
shows that {idzQa,b, i^) + M{i^ + y) = 0. Thus = i'iQa,b + i^2^dzQa,b-, where 

{Qa,b,Qa,b) {Qa,b,idzQa,b) \ ( _ _( ^i^ + y) 

(}dzQa,b,Qa,b) {^dzQa,b,idzQa,b) J \V2 ) \M{v + y) 
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Observe that the matrix of this system is exactly the same one as in ()3.17|) . Thus, proceeding as 
above, we obtain \\y\\x < {C/a)\\i' + some C > independent of a,h. Since we aheady 

know that \\v\\x 1^ Cina, it fohows that \\i'\\x 1^ {C2 / a)\\y\\'j^ provided k > is sufficiently small. 
□ 

To show that the energy £a,b has a local minimum on S^fc at Qa,bj we consider the second 
variation of Sa,b at Qa,b, i-e. the linear operator Ha^b defined in H3.13() . 

Lemma 3.7 // ||(a, is sufficiently small and ah 7^ 0, then 



Proof. We use a perturbation argument. When (a, = (0,0), the operator Ha^b reduces to a 
differential operator with constant coefficients: Hq = —Adzz — ^^dz- This operator is self-adjoint 
in the real Hilbert space Lpgj,([0, 27r], C) equipped with the scalar product (•,•), and its spectrum 
is (t{Hq) = {4n(n it 1) | n G Z}. The null space of Hq is spanned by the four vectors ^0, i^Oi 
where ^0 = and r]o = 1 (see ()3.15|) ). The other eigenvalues of Hq are positive and greater or 
equal to 8, hence the quadratic form ho : X ^M. associated to Hq satisfies 



We now consider the quadratic form ha^b ■ X ^ defined by ha^biQ) = {Ha,bQ,Q)- This form 
is uniformly bounded for (a, b) in a neighborhood of zero, i.e. there exists C3 > such that 
ha,biQ) < CsWQfx for all QgX. Moreover, h 

a^b converges to /iq as (a, b) — > (0, 0) in the following 

sense: 



In particular, we have ha^biQ) ^ ^HQH^ for all Q £ Yq if ||(a,6)|| is sufficiently small. 

On the other hand, since \\Ca,b — ^o\\x + WVafi — "HoWx = 0{a'^ + 6^), it is straightforward to verify 
that the subspace 



{Ha,by,y) > 6||y||x , for all y G . 



(3.18) 



ho{Q) = {HoQ,Q) > 8\\Q\fx , for ah QeYo 



where 





Ya,b = {Q &X\{Ca,b,Q) = (iCa,fe,Q) = {r]a,b,Q) = {^r]a,b,Q) = 



converges to Yq as (a, b) — > (0, 0) in the following sense: 



S{Ya,b,Yo) =^ sup{distx(Q,ro) I Q G Ya,b , WQWx = 1} = 0{a^ + b^) . 



In particular, if Q G Ya^b satisfies HQHx = 1, we can find Q G Yq with ||Q||x = 1 and \\Q — Q\\x as 
small as we want, provided (a, b) is close to zero. Since hafi{Q) > 7 and 



ha,b{Q) — ha,b 



iQ)\ < \\ha,b 



\\{\\Q\\x + \\Q\\x)\\Q-Q\\x < 2C3\\Q-Q\\x 



we conclude that ha^b{Q) > 6 if ||(a, &)|| is sufficiently small. This proves ()3.18() . 



□ 



Using Lemmas 13.51 and 13. 7( we are able to give a more precise version of Claim 1 above: 
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Lemma 3.8 There exists ki > such that, if ||(a, b)\\ is sufficiently small and b > a > 0, then for 
all Q G Tia,b satisfying p{Q,Qa,b) < i^io- one has the inequality 

£a,b{Q) - £a,b{Qa,b) > p{Q , Qa,bf ■ (3.19) 

Proof. If Q € '^a,b satisfies p{Q,Qa,b) ^ for some ki > sufficiently small, we know from 
LemmaEinithat Q = 1l(^,^)iQa,b+i^+y), where {ip,(,) £G,u £ Na,b, y G i^a.b, \\y\\x < Cip{Q,Qa,b), 
and < (C'2/a)||y|lx- particular, < «^iC'iC2||y||js:- Since the energy Sa^b is invariant 

under the action of G, we have £a,b{Q) = £a,b{Qa,b + v + y). As £'^j,{Qa,b) = and £'^f^{Qa,b) = Ha,b, 
we obtain using Taylor's formula: 

£a,b{Q) - £a,b{Qa,b) = ^iHafiiy + 1^), {y + 1^)> + 0{\\yfx) ■ 

But {Ha^by,y) ^ 6||y||^ by Lemma 1^71 hence 

^{Ha,b{y + ^),{y + ^)) > 3||y||x - C'sllylUllz^llx - ^C'slli/ll^ 

> l|y||x(3 - KiCiCsCs - U'^iCiC2fC3) , 

where C3 is the constant in the proof of Lemma 13.71 Thus, if ki is sufficiently small, we obtain 
£a,b{Q) — £a,biQa,b) > Under the same assumption, we also have p{Q,Qa,b) < \\y + < 

\\y\\x{l + K1C1C2) < V2\\y\\x, and (1X1111) follows. □ 

3.4 Proof of Proposition 13.21 

The proof of Proposition 13.21 consists of three steps in which we treat successively the three types 
of waves: the zero solution (a = b = 0), the plane waves {ab = 0), and the periodic waves {ab 7^ 0). 
In each case, the arguments rely upon energy estimates as the one given in Lemma 13.81 for the 
periodic waves {ab 7^ 0). In the case of the plane wave (5o,fe = b the stability is proved using the 
functional 

em = m) - b'N{Q) = ((2-362)|Q,|2 + _(|Q|2 _ ft2^2j dz - — , 
for which we have the analog of Lemma 13.81 

Lemma 3.9 There exists K2 > such that, if b > is sufficiently small, then for all Q G X 
satisfying p{Q, Qo,b) ^ 1^2^ and N{Q) = A^((3o,fe) one has the inequality 

SbiQ) - £b{Qo,b) > \p{Q,Qo,b? ■ (3.20) 
o 

Proof. Without loss of generality, we assume that ||Q — Qo,fe||x < ^2^. Since 

\Q{z) -b\< \\Q - Qo.felU <K2b<b , for all z G M , 
we can write Q = {b + r)e^^ , where r, 99 S X are real functions satisfying 

\r{z)\<n2h, |e''^('^) - 1| < 2k2 , for aU z G M . 
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These inequalities imply ||r||x + ftllv^Hx < CK26, for some C > independent of a,b. Thus 
SbiQ) - £biQo,b) = l^^{{2-3b'){rl + {b + rf^l)+^{2br + r^f]dz 



> r(rl + b\l)dz. (3.21) 







On the other hand, since N{Q) = iV((5o,fe)) "we have J^^ {2br + r^) dz = 0, and Poincare's inequality 
implies 

I 1 + TTT I dz < / ( 1 + - I dz , hence / dz < 2 rldz . 



^ 26/ Jo -/o -/o 

Finally, \i (p = {2-k)~'^ J^^ vi^) dz, we have 

p(Q,Qo,b)' < ||Q-6e''?||^ < 2||re''^||2, + 262||e''^-e'^||^ 

/'27r /'27r 

< 2||r||^ + 2/ rV^ dz + 26^1199- (^11^ < 6 / (r^ + 6V^)dz, (3.22) 

JO 

again by Poincare's inequality. Combining (|3.21j) and (|3.22j) . we obtain (|3.20j) . □ 

We are now in position to prove Proposition 13.21 

Proof of Proposition 1231 Throughout the proof, we assume that ||(a,6)|| is sufficiently small 
and that 6 > a > 0. Given Qq ^ X with p{Qo,Qa,b) < we consider the solution Q{z,t) of 1)3. 4() 
with initial data Qq. Replacing Qq with 7^(^^^)Qo if needed, we can assume that ||<5o — Qa,6||x < £• 
We distinguish three cases: 

Case 1: a = 6 = 0, i.e. Qa.fe = 0. In this case, if e > is small enough, the solution Q{-,t) of l|3.4|) 
satisfies ||(5(-,t)||x < 2e for all t G M. This is obvious because the quantity 

£{Q)+m{Q) = ^'"(2|g,|2 + i|Q|4 + 2|Q|2)dz 

is invariant under the evolution of ((231), and satisfies 2||Q||^ < £{Q) + 4:N{Q) < 4||Q||^ if ||Q||x 
is small. 

Remark: As a consequence of this preliminary case, we assume from now on that e < K3(a2 + 62)1/2 
for some small K3 > 0. Indeed, if e > ^3(0^ + 62)1/2^ ^j^g trivial estimate 

\\Q{;t)-Qa,b\\x < \\Q{;t)\\x + \\Qa,b\\x < 2\\Qo\\x + \\Qa,b\\x <2e + ^Qa,b\\x , 
which gives the desired result since HQa.fcHx < C{a^ + b"^)^/^ < {C / K^)e. 

Case 2: 6 > a = 0, i.e. Qa,b = 6 is a plane wave. We consider initial data Qo ^ X such that 
IIQo - QoAx < e < H3b. If iV(Qo) = N{Qo,b), then N{Q{-,t)) = N{Qo,b) for all t G M, and 
Lemma 13.91 implies 

p{Q{;t),Qo,bf < (^{£b{Q{;t))-£b{Qo,b)) = Q{£b{Qo) - £b{Qo,b)) < Ce^ , 
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provided that Ce^ < KqP'^, which is the case if Cn^ < k^. If N{Qq) ^ A^((3o,fe)) we define m = 
TT~'^iN{Qo)- N{Qo^b)), so that N{Qo) = iV(Qo,fe)! where Q'^,, = (b^ + uj)^^^. So we are led to study 
the stabihty of the rotating wave Qg^^e"''^* of (|3.4() with respect to perturbations preserving the 
charge A^. This can be proved exactly as above, and we obtain p{Q{-,t),Q^^) < Ce for all t G M. 
Since \\Qqi, — Qo,b\\x — C\u!\/b < Ce, we have the desired result. 

Remark: As a consequence, we can assume from now on that e < n^a for some small K4 > 0. 
Indeed, if e > K^a, we can use the easy estimate 

p{Q{-,t),Qa,b) < piQ{-,t), Qo,b) + \\Qo,b - Qa,b\\x ■ (3.23) 

Observe that \\Qo,b — Qa,b\\x < Ca for some C > independent of a, b. In particular 

IIQo - Qo,b\\x < WQo - Qa,b\\x + \\Qa,b " Qo,b\\x < e + Ca < {1 + C/K4)e , 

hence for e > small enough we have p{Q{-,t), Qo^b) ^ Ce for all t G M. It then follows from H3.23|) 
that p{Q{--,t),Qa,b) < C"e for all t G R, which is the desired result. 

Case 3: 6 > a > 0, i.e. Qa,b is a nontrivial periodic equilibrium of 1)3. 4j) corresponding to a 
periodic wave of (|1.2|) . Assume that Qq ^ X satisfies \\Qq — QafiWx < e < h^cl- If <5o £ ^a,6) then 
Q{-,t) £ ^a,b for alH G M and Lemma 13.81 implies that 

p{Q{-,t),Qa,b? < £a,b{Q{-,t))-£a,b{Qa,b) = £a,b{Qo) - £a,b{Qa,b) < Ce^ , 

provided Ce^ < k^o?, which is the case if Ck\ < k\. If Qo ^ '^a.b, then by Claim 2 above there 
exists (a;,c) G with |a;| + |c| < Cbe such that N{Qq) = NiQ'^'f) and M(Qo) = Af(Q^p. So we 
are led to study the stability of the periodic wave u{x, t) = QHPa,bX^i-i-+^)t) Q'^'^(2ka,bX + ct) of (tT^ 
among solutions of the form e^^P'^-''^~^^ Q{2ka^bX,t) for which the charge N and the momentum 
M have the same values as for the periodic wave. But if we apply a dilation of factor A and 
a Galilean transformation of speed v, where X,v are given by ()3.9() . the periodic wave becomes 
u{x,t) = e^^^°'''''^~*^ Qa'fi'C^ka'^b'x) for some {a',b') close to (a, 6), and we are back to the previous 
case. As \\Qo-Qa',b\\x < WQo - Qa,b\\x + \\Qa,b - Q^'^bWx < Ce, this shows that t), Q^p < Ce 
for all t G M, and the result follows. This concludes the proof of Proposition 13.21 □ 

Remark 3.10 In \15f the authors use the decomposition X = Ta^b ®-^a,b, where N'a,b is the two- 
dimensional space spanned by 



d 



d 



duiQafi — Q^Qa'b r ...^_rnr.^ ' ^cQa,b — r:^^Q a]b 



{w,c)={0,0) OC 



(c^,c)={0,0) 



This alternative decomposition has the advantage that {Ha^b^, = for all u G Ma,b Oi'^d all 
V G Ta,b, because 

Ha,b{d^Qa,b) = N'{Qa,b) = Qa,b , Ha,bidcQa,b) = M'{Qa,b) = id,Qa,b ■ 

Using in addition i'J.ll]) we also find 

^ _ _f {Ha,b{du>Qa,b), dujQa,b) {Ha,b{dujQ a,b) dcQ a,b) 
\ {Ha,bidcQa,b),dujQa,b) {Ha,b{dcQ a,b) , ^cQ a,b) 
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Remark that the spaces Ma,b cL^d Ma,b cli^g very close when (a, b) is small, since 

du.Qa,b = 7^(a + be'''){l + 0{a^ + b^)) , dcQa,b = ^{2a - 6e-'^)(l + 0{a^ + b^)) . 
bab bab 

Remark 3.11 (2n7r— periodic perturbations) As an intermediate step between the periodic set- 
up considered in Theorem Q and the case of arbitrary bounded perturbations for which no result is 
available so far, one can try to study the orbital stability of the travelling waves of il.^) with respect 
to perturbations whose periods are integer multiples of the period of the original wave. This amounts 
to replacing the space X in Provosition VJ.'/^ by H^^^{\}), 2n7r], C) for some n > 2. In that case most 
of the results above remain valid, but the linear operator Ha^b has now 2n — 1 negative eigenvalues. 
Thus the number of negative eigenvalues of Ha^b minus the number of positive eigenvalues of the 
Hessian matrix 'Ha,b is equal to 2n — 2, a nonzero even integer. This means that neither the Stability 
Theorem nor the Instability Theorem in J15f applies if n > 2. The only way out of this difficulty 
would be to replace the manifolds Tia^b defined in ( |5'.7| ) by invariant manifolds higher codimension, 
which amounts to use additional conserved quantities of instead of N and M only. 



4 Spectral stability 

In this section we prove the spectral stabihty result in Theorem |2j We start with the evolu- 
tion equation ()3.4() . which we linearize about the stationary solution Qa,b{z) corresponding to the 
periodic wave Ua,b{x,t) = e^^^'^'''^~^^ Qa,bi2ka,bx) of the NLS equation. We find the linear operator 

Aa,bQ = ^^klbQzz - ^Pa,bka,bQz + i(l " plb)Q " 2i|Qa,6pQ " '^Q^bQ , (4.1) 

which we consider in either the real Hilbert space Y = L^(M, C) (localized perturbations) or the real 
Banach space Y = Cb(M,C) (bounded pertubations). To study the spectrum Aa,b^ it is convenient 
to decompose the elements of Y into real and imaginary parts, in which case we obtain the matrix 
operator 

, ^( -^Pa,bka,bdz + 2Ra,bIa,b 'KA^ + Kb " 1) + ^b + Kh 

V ^^lA- - Kb - 1) - ^Kb - ^a,b -^Pa^bkaA " 2Ra,bIa,b 

where Qa,b = Ra,b + ^Ia,b- We are now interested in the spectrum of this matrix operator in the 
(complexified) spaces L^(M, C^) and Cb(R, C^). We prove that the spectrum of Aa,b in both spaces 
lies entirely on the imaginary axis, if ||(a, 6)|| is sufficiently small. This means that the periodic 
wave C/a^b is spectrally stable in Y . 

4.1 Bloch-wave decomposition and symmetries 

The spectral analysis of Aa,b relies upon the so-called Bloch-wave decomposition for differential 
operators with periodic coefficients. This method allows to show that the spectrum of Aa,b is 
exactly the same in both spaces L^(M, C^) and Cb(M, C^), and can be described as the union of 
the point spectra of a family of operators with compact resolvent (see e.g. |2U[ I22j ). In our case, 
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the operator Aa,b lias 27r-periodic coefficients and its spectrum in both L^(M, C^) and Cfe(M, C^) is 
given by 

(T{Aa,b) = (j f^(A,fe,7) ' (4-3) 

where the Bloch operators 

, ^( -"^PafikaAdz + i7) + 2Ra,bIa,b -^b^^z + h? + (p^fe " 1) + Rib + 

V^^'-''^^^ + i^)' - ^Plb - 1) - 3Rlb - lib -^Pa,bka,b{d. + i7) - 2R,,bIa,b ) 

are hnear operators in the Hilbert space of 27r-periodic functions Lpgj,([0, 27r], C^). We can now 
reformulate the spectral result of Theorem |51 as follows: 

Proposition 4.1 There exists (5i > such that, for any 7 G (~^) |] o.'^T'd any {a,b) € with 
||(a,6)|| < 61, the spectrum of the operator Aa,b,'y in Lpgj.([0, 27r], C^) satisfies (7{Aa,br() 



We equip the Hilbert space Lp(,j.([0, 27r], C ) with the usual scalar product defined through 
((Qi, g2)*, {Ri,R2Y) = j^^ (Qi{z)Ri{z) + Q2{z)R2{z)) dz . 



10 

The operators ^a,6,7 are closed in this space with compactly embedded domain //pgj.([0, 27r], C^). 
An immediate consequence of the latter property is that these operators have compact resolvent, 
so that their spectra are purely point spectra consisting of isolated eigenvalues with finite algebraic 
multiplicities. Our problem consists in locating these eigenvalues. 

The spectra of the operators Aa,b and Aa,b,-y possess several symmetries originating from the 
discrete symmetries of H1.2() and the symmetries of the wave profile Qa,b- First, since the operator 
Aafi has real coefficients, its spectrum is symmetric with respect to the real axis: cr{Aa,b) = (^{•^a,b)- 
For the Bloch operator Aa,b,'y, the corresponding property is cr{Aa,b,'r) = cr{Aa,b,-'y)- Next, it is 
straightforward to check that Aa,b has a reversibility symmetry, i.e. it anticommutes with the 
isometry S defined by 

Qi{z)\ [ Qi{-z) 



Q2{Z)) \-Q2{-z)J ■ ^^-^^ 

Thus SAa,b = —•^a,b<S, which implies that the spectrum of Aa,b is symmetric with respect to 
the origin in the complex plane: cr{Aa,b) = —c{Aa,b)- The corresponding property for the Bloch 
operators is SAa,b,-y = —Aa,b-'y<S, which implies that cr{Aa,b,'y) = —o'{Aa,b,-'r)- Summarizing, the 
spectrum of Aa,b is symmetric with respect to both the real and the imaginary axis, and the spectra 
of the Bloch operators Aa,b,'y satisfy 



In particular, the spectrum of .4a, 6,7 is symmetric with respect to the imaginary axis and we can 
restrict ourselves to positive values 7 S [0, ^] without loss of generality. 

Using now the relations (|3.2|) for the wave profile Qa,b, we see that the spectra of Aa,b satisfy 
(^{•^a,b) = (^{•^-a,b) = '^{•^-a-b) and o-{Ab^a) = —f^i-^afi)- (Actually, the last equality is easier to 
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establish if we use the complex form (|4.1I) of the operator Aa,b, for which we have Ab^ai^ ^^Q) = 
—e~^^Aa,bQ-) Similarly, we find for the Bloch operators 

(^{Aafin) = (^i-^-aAl) = ^(-^-0,-6,7) , and a{Ab,a,-y) = -(^{Aa,b,-y) ■ (4.6) 

Finally, we note the formal relation Aa,b = —^Ha^b between the linearized operator ()4.1|) and 
the second variation of the energy defined in H3.13|) . For the Bloch operators, we can write in a 
similar way 

with 

^ l-K,bidz + + iplk - 1) + ^K,b + 4% ^Pa,bka,b{d. + i7) + 2Ra,bIa,b \ 

V -'^PaMdz + i7) + 2Ra,bIa,b "4/^2 ,(9, + i7)2 + (pI, - 1) + Rl, + 3/2 , J ■ 

Actually, this property is a consequence of the Hamiltonian structure of the NLS equation. Though 
some properties induced by this structure are exploited, we shall not make an explicit use of the 
Hamiltonian structure itself in the proof of spectral stability. 

4.2 First perturbation argument and properties of the unperturbed operators 

Our spectral analysis for the operators ^a,fe,7 relies upon perturbation arguments in which we 
regard Aa,b,-y as small bounded perturbations of the operators with constant coefficients 

,0 ^ ( -'^Pa,bka,b{dz + i7) -4A;2 + i7)2 + {pi, - I) \ 

The difference A\f^ := ^a,6,7 — -^067 is a bounded operator with norm H^^^H = 0{a? + 6^), as 
(a,6)^(0,0). 

A straightforward Fourier analysis allows to compute the spectra of the operators ^4^^^^^: 

^^(^^,7) = {^^ib,v 4ha = -^PaMn + 7) ± (4A:^,,(n + 7)^ + p^^, - l) , n G Z} C iM , (4.7) 
in which the eigenvalues are all semi-simple with eigenfunctions 

Furthermore, the resolvent operators '^^^-^(A) = (Al — A^^f^^)~^ have norms 

ll<^7(^)ll = dist(A,0-K,,7)) ' ^ ^ "'^^ ■ 

A simple perturbation argument shows now that the spectrum of ^a,6,7 stays close to o"(^^^^) 
provided ||(a, is sufficiently small. More precisely, we have the following result. 
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Lemma 4.2 For any c > there exists 6 > such that for any 7 G [0, ^] and any (a, b) G with 
II (a, 6) II < 5 the spectrum of Aa,b,-y satisfies 

a{Aa,b,^) C IJ B{iu;-f^^;c) U \J B(i<^^;c), 
in which B{iuj^'^^;c) represents the open ball centered at iw^^"^ with radius c. 

Proof. For any A ^ U„6Z ^(i^^Cy' 

Al - = - ^° fe,^ - = (Al - A%^^) (1 - n%^^{X)Al,) . 

Since 

upon choosing 6 sufficiently small, we have that 1 — TZ^ jj^{X)A\ is invertible, so that Al — Aa,b,'y 
is invertible, as well. This proves that A does not belong to o-{Aa,b,'y)- 

In order to locate the spectra of Aa,b,'r, we need a more precise description of the spectra of 
A^ J ^. Looking at a = 6 = we find that: 

• if 7 = 0, all nonzero eigenvalues of q q are double, 

+,n +,1— n — ,n — , — 1— ra 

'^0,0,0 ~ ^^0,0,0 ' ^^0,0,0 ~ ^^0,0,0 ' 



and zero is an eigenvalue of multiplicity 4, 



±,0 _ +,1 _ -,-1 _ p, 
^0,0,0 — ^^0,0,0 — ^^0,0,0 — u 



if < 7 < ^ all eigenvalues are simple; 

if 7 = ^, there is a pair of simple eigenvalues iti. 



--1 +,0 1 

^0,0,1 = -^0,0,1 = 1 ' 

and the other eigenvalues are all double, 

+,n +,— n — ,n —,—2—n 

0,0,i 0,0,i ' 0,0,i 0,0,i 

We therefore distinguish three cases: 7 Ri 0, 7 |, and 7 e [7*, | — 7*] for some 7=,, G (0, |), 
which we treat separately in the next paragraphs. In each case, the starting point is an estimate of 
the distance between any pair of eigenvalues of A^ which is directly obtained from the explicit 
formulas H4.7|) . We use this estimate to construct an infinite family of mutually disjoint sets (balls 
or finite unions of balls) with the property that the spectrum of Aa,b,'y is contained in their union. 
Inside each set Aa,b,'y will have a finite number of eigenvalues (one, two or four) so that the problem 
reduces to showing that these eigenvalues are purely imaginary. In Propositions 14.41 14. 6( 14.81 and 
I4.1UI below, we show that in all three cases the spectrum of Aa,b,'y is purely imaginary, provided 
II (a, 6)11 is sufficiently small. This proves Proposition 14.11 
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4.3 Spectrum for 7 away from and ^ 

We start with the case 7 G [7*, ^ — 7*], when the operators with constant coefficients -4.^^^^ have 
only simple eigenvalues: 

Lemma 4.3 For any 7* G (0, j), there exist positive constants and such that for any 7 G 
[7*,^ —7*] and any {a,b) with ||(a,6)|| < 5*, we have 

\^^a'b'y ~ ^^at-yl — ^* ' f'^'^ all p,n G Z and all a,T e { — ,+} with {a,n) / {t,p) . 

This lemma shows that the distance between any pair of eigenvalues of b is strictly positive, 
uniformly for small ||(a,6)|| and 7 G [7*,^ — 7*]. This allows us to find an infinite sequence of 
mutually disjoint balls with the property that the spectrum of Aa,b,'y is contained in their union, 
and that inside each ball both operators have precisely one simple eigenvalue. The symmetry of 
the spectrum of ^a,fe,7 with respect to the imaginary axis then implies that this simple eigenvalue is 
purely imaginary. We point out that classical perturbation results for families of simple eigenvalues 
jl8j do not directly apply, since here we have infinitely many eigenvalues. 

Proposition 4.4 Fix 7* G (0,^). Then there exist positive constants a and 6 such that for any 
7 G [7*, I — 7*] and any (a, b) with ||(a, 6)|| < 6, the following properties hold: 

(i) The spectrum of Aa,b,'y satisfies 

^^(Aat) C U B{iu^:M,c) U U i3KX;c) , 

and the closed balls S(iu;^^"^;c) are mutually disjoints. 

(a) Inside each ball B{iuj^'^^;c) the operator Aa,b,'y h-as precisely one eigenvalue, which is purely 
imaginary.^ 

Proof. (i) We may choose any c < c*/4 with c=k the constant in Lemma 14.31 so that the balls 
B{iu>^'i^^; c) are mutually disjoints, and then apply Lemma 14.21 

(ii) Consider a ball i?(ia;^^"^; c). Inside this ball has precisely one eigenvalue ioj^'^^ with 

associated spectral projection n[|'^^ satisfying ||n||'^^|| = 1. The result (i) provides us with a 
spectral decomposition for the operator Aa,b,'y^ and we can compute the spectral projection 11^^^^ 
associated to B(iu^'^^^;c) via the Dunford integral formula 

Km = ^a,6,7(A)dA , (4.8) 

in which C„ is the boundary of B{\oj^'^^; c) and TZa^b,'yW = (^^1 — •^0,6,7)"^- Using the formula for 
the resolvent 

^^^^^^(A^<,,^(A) (1 - <,<,,^(A))-' , 

^Here and in the rest of the paper we say that "yla,^,^ has n eigenvalues inside the set B" when the sum of the 
algebraic multiplicities of the eigenvalues of Aa,b,-y inside B is equal to n. 
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which holds for sufficiently small 6 since 11-4.^6 1| = 0(a^ + ^^), we compute the difference 



Since ||'^afe^(A)|| = 1/c, for A G Cn, we conclude that 

lin:..-n-|,.E(l|K.|,)^^ 



Upon choosing 5 small enough we achieve 



|n^.,-n^'" Jl < 



0,n 1 



a,b,7 a,b,-y\\ 1 i ||n" _ n°'" II IITT^'" II _l_ lln" _ n^'" I 



< min 



1 



IttO," II ' IITT" I 
l-^-^a,6,7ll II a.,b,'r\ 



so that the projections HJ^j^^ and n|^'^^ realize isomorphisms between the associated spectral 
subspaces of A^f^^ and Aa,b,'Y (1^1 Lemma B.l]; see also ^1 Chapter I §6.8]). In particular, 
they have the same finite rank, so that Aa,b,'y has precisely one simple eigenvalue inside the ball 
B{iLJ^'^^; c). Finally, since the spectrum is symmetric with respect to the imaginary axis (|4.5() this 
simple eigenvalue is necessarily purely imaginary, which concludes the proof. □ 

4.4 Spectrum for small 7 

We start again by analyzing the distance between the eigenvalues of A^^^^, now for small values 
of 7. Since at = 6 = 7 = the spectrum of A'lf^^ consists of double nonzero eigenvalues and a 
quadruple eigenvalue at zero, for small a, b, and 7 we expect pairs of arbitrarily close eigenvalues 
together with four eigenvalues close to the origin. A precise description of the location of these 
eigenvalues is given in the following lemma. 

Lemma 4.5 There exist positive constants 70, c, and 6, such that the following holds, for any 
7 G [0, 7o], and {a,b) with ||(a, < 6: 

(^) i^S7' '"^ttv ^"^aM ^ B{0;1); 
(U) iu;-',\, ia.+-\ iu^ i B{0; 4), n G Z \ {-1, 0, 1}; 
(Hi) - ^^afj ^ c, G Z \ {0, 1}, p / p / 1 - n; 

(^^) li'^o^Ci ~ ^^oiCi' >c,n,peZ\ {-1,0}, p^n, -1-n; 

M - >c,nez\{0,i}, pez\ {-i, 0} . 
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The first two properties (i)-(ii) in this lemma together with the perturbation result in Lemma [4.2l 
provides us with a spectral splitting for Aa,b,'y- 

) = 0'i(^a,(),7) ^ (-^0,6,7) , 

with 

(•^a,b,7 

) C 5(0; 2) , (72(^,6,7)^5(0; 3) = . 

Inside the ball 5(0; 2) we find the part of the spectrum of Aa,b,'y which is close to the quadruple 
zero eigenvalue of ^qqq, whereas the rest of the spectrum lies outside the ball 5(0; 3). The last 
properties (iii)-(v) show that the eigenvalues outside 5(0; 3) are well separated except for the pairs 

(i'^^C7'^^'t^7")' ^ ^ ^ \ {0' 1}' and (i^^";"^, it^~;77^~"), ^ ^'^\ {-I'O}' which may be arbitrarily 
close. At a = 6 = 7 = 0, these are precisely the double eigenvalues of Aq q q- Notice however that, 

for fixed a, b, 7, the distances |iw^^"^ — ^^ab'y'^\ ^^'^ \^^ab''y ~ ^^a'b~-^~^\ typically grow like 0{'n?) 
as |n| ^ cxD, see 1)4. 7() . 

We analyze these two parts (71(^0,6,7) and cr2{Aa,br() of the spectrum of ^a,6,7 separately in the 
Propositions 14.61 and 14.81 below. 

Proposition 4.6 There exist positive constants 70, c, and 5, such that for any 7 G [0, 70], and 
{a,b) with II (a, 6) II < 6, the following holds. 

(i) The spectrum a-2{Aa,b,'r) satisfies 

a2(A,6,7)c5(iu;-;^;c)U5(i^+-;;c)U [j 5(iu;-;;c)U [j B{iu;X;c) , 

n^±l,0 nf^±l,0 



in which the balls B{iLO^'^^;c) are mutually disjoints, except for some pairs (iw^J^^, iw^^^^ 
neZ\{0,l}, or (iu;-^,'i^--i-"), nEZ\ {-1,0}. 



(ii) Inside each ball B{iuj^'^^;c) the operator Aa^b,-y has either one or two eigenvalues, which are 
purely imaginary. 

Proof. The result (i) is obtained from Lemma 14.51 (iiiVfvl and Lemma 14.21 just as the first part 
of Proposition 14.41 The only difference is that here we have pairs of balls which are not disjoint. 
As was noticed above, the distances |iw^^"^ — i'^^fe^-^^l and |iw~^"^ — i<^afe~7^~"l grow like 0{n'^) as 
|n| — > 00, so that we have in general a finite number of such pairs for a given value of a, 6, 7. 

(ii) For the balls 5(ia;^^"^;c) which are disjoint from all the others we can argue and conclude 
as in the proof of Proposition 14.41 It remains to consider the case of two balls which are not 
disjoint. Choose a pair of eigenvalues (io;^^"^, ia;^^^~") such that B(iijj~^'^^;c) n B{\lo'^'^~'^;c) 7^ 
(the argument is similar for a pair {loJ~'^^,\uJ~'^^^~"')). We construct the spectral projection n"'^~" 
for Aa,b,'y corresponding to the union of these balls with the help of the Dunford integral formula 
(|4.8|) . in which the circle Cn is replaced by the smallest circle Cr with radius c < r < 2c, centered 
on the imaginary axis, which contains both balls. The spectral projection IlJJ'?''^"" for A^^,^ has 
unit norm again, and since H'/^a ^(A)|| < 1/c for A € Cr, we easily find 



\\Ai II 2\\A'^ I 

|n.M - n,,,, II < - ^^prj ^ ^^p^ = 0(a + 6 ) . (4.9) 
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As in the proof of Proposition 14.41 we now choose 5 sufficiently small such that these projections 
realize isomorphisms between the associated spectral subspaces oi A^^^^^ and Aafi^^y- In particular, 
these subspaces have the same finite rank equal to 2, which proves that Aa,b,'y h^-s precisely two 
eigenvalues in B{\u'^^^^^] c) U 5(iw^fe^~"; c). 

In order to show that these two eigenvalues do not move off the imaginary axis we choose 
an appropriate basis of the associated two-dimensional eigenspace and compute the 2x2 matrix 
representing the action of Aa,b,'y on this space. Then it suffices to show that this matrix has purely 
imaginary eigenvalues. We start with the basis 

of the two-dimensional eigenspace of -4.1^^7' which satisfies {—'^J£,Qk-,im) = ^M- We claim that for 
-^a,fe,7 we can find a basis with the same property. Indeed, consider the vectors 

which form a basis of the two-dimensional eigenspace of Aa,b,'y From (|4.9j> we obtain (— iJ^Oi '^o) = 
1 + 0(0^ + 6^) > 0, so that the vector ^0 defined by 

^° " / • j7 7x112 ^° ' 



satisfies (— iJ^OiCo) = 1- Then we define successively 



and find (— i^CoiCi) = and (— iJ^iiCi) = li which proves the claim. 

The property {—iJS,k,(,£) = implies that the action of ^a,b,7 on the two-dimensional space 
spanned by {^016} is given by the matrix 



M 



a, 6, 7 



Using the decomposition Aa^b^^y = JHa,b,'r we find 

{Aa,b,'yCk,-^JW = {'JHa,b,-iik-,-''^Jil) = {Ha,b,'yS,k, -^J~^ J^l) = i{Ha,b,'yCk , , 

SO that 

_ f {Ha,b,'rCo,Co) {Ha,b,'rCl,Co) 



Since {Ha^b,-yQ, R) = {Ha,b,jR,Q), we conclude that this matrix always has purely imaginary 
eigenvalues. This completes the proof. □ 



Remark 4.7 The last part of this proof is a simple version of the well-known result for general 
Hamiltonian systems which asserts that colliding purely imaginary eigenvalues do not leave the 
imaginary axis when they have the same Krein signature (see e.g. |i7| /j. In the case of the four 
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eigenvalues close to the origin, which we treat in the next proposition, the same argument does not 
work anymore (these eigenvalues have opposite Krein signatures) . Instead, we compute an explicit 
expansion of the restriction of ^a,fe,7 the associated eigenspace which allows to show that these 
four eigenvalues are purely imaginary. 

Proposition 4.8 There exist positive constants 70, c, and 5, such that for any 7 G [0, 70], and 
{a,b) with ||(a,6)|| < 6, the set ai{Aafi,"/) consists of four purely imaginary eigenvalues. 

Proof. As in the previous cases, upon choosing 6 sufficiently small, we obtain that Aafi,"/ 
has precisely four eigenvalues inside the ball -6(0; 2). In order to locate these four eigenvalues 
we construct a suitable basis for the associated eigenspace and compute the 4x4 matrix M.a,b,'y 
representing the action of Aa,b,'r on this space. Then we show that this matrix has purely imaginary 
eigenvalues. 

We start with the particular cases a = b = 0, and 7 = 0. In the first case, the operator 
-^0,0,7 = -^0,0,7 constant coefficients, so that we can explicitly compute the basis and the 
matrix. We choose the real basis 

(0) _/sinz\ (1) _ /0\ (2) _ / COSZ \ (3) _ /I 

^O'O'^ " \cosz) ' ^O'O'^ " Vly ' ~ \-smz) ' ^O'"-^ " Vo 

in which we find 

Next, we consider the operator Aa,b,o- As for the operator Ha^b in Section 13 we have that 



d_ 



d 

T^{ip,^)Qa,b = —^Qa,b , -^'^(^,^)Qa,b 



, , ^ ~ 9zQa,b , 



(^,0=(o,o) ' 
belong to the kernel of Aa,b,o- In addition, since 

Ha,b{d^Qa,b) = N'{Qa,b) = Qa,b , H,,b{dcQa,b) = M'{Qa,b) = id.Qa^b , 

(see Remark l3.1U() and Aa,b,o = —^Ha^b, we have 

Aa,b,o{dujQa,b) = —^Qa,b , Aa,b,o{dcQa,b) = dzQa,b , 

which provides us with two principal vectors in the generalized kernel of Aa,b,o- Together with the 
two vectors in the kernel of Aa,b,o these give us a basis for the four-dimensional eigenspace of Aa,b,o- 
At a = 6 = we must find the basis above so that we set 

= l(.iQa,b + dzQa,b) = (°) +0(|a|(|a| + |6|)) , 
o = 2ai2d^Qa,b - dcQa,b) = (.""g-^^J +0{a^ + b^) , 
^AO = '^Kdu.Qa,b + dcQa,b) = +0(a2 + 62) , 
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and a straightforward calculation gives the matrix 

/O9 M2(a,6)\ /O 0\ /-2o2 -4a6\ ,^,4 

Finally, we consider the full operator Aa,b,-y and construct a basis {C^al'y^^a^l'y^^'ab'y^^'ab'y} 
small a, 6, and 7, by extending the bases above. Notice first that the vectors in the basis for 7 = 
satisfy 

- - 9^(2) _ t(2) Cd(3) _ ^(3) 

where 5 is the reversibility operator H4.4I) . Since for 7^0 we have SAa,b,'y — ~-^a,b,—'f'S , the 
vectors in the basis can be taken such that 

"^0,6,7 ?a,6,— 7' "^0,6,7 "^0,6,— 7' ^a,b,7 ?a,6,— 7' ?a,6,7 ?a,b,— 7 ' 



and then the matrix A4a,b,"f satisfies 

SMa,b,'y = -Mafi-'yS , where S 



-I2 O2 
O2 I2 



In addition, since Aafi = A-a-b-, we also have A1a,fe,7 = -^-0,-6,7- Together with the results for 
a = 6 = and 7 = we conclude that 



Ma,ba 



4i7 (D2 + 0{a^ + 62)) M2(a, h) - 472 (I2 + 0{a^ + ^2)) 
472 (I2 + 0(a2 + 62)) 4i7 (D2 + 0(a2 + 62)) 



To end the proof we show that the four eigenvalues of this matrix are purely imaginary. The 
structure of the matrix A^a,fe,7 implies that its characteristic polynomial is of the form 

+ i7C3A^ + 7^C2A^ + i7^ciA + 007^^ , 

in which the coefficients Cj depend upon a, 6, and 7. The four roots of this polynomial are 
symmetric with respect to the imaginary axis, because the spectrum of ^a,fe,7 is symmetric with 
respect to the imaginary axis, so that the coefficients Cj are real functions of a, 6, 7. In addition, 
the spectral equalities ()4.5|) " H4.6p imply that Cj are even in a, 6, and 7, and that when replacing 
(a, 6) by (6, a) the coefficients cq, C2 do not change, while ci, C3 change sign. 

We now set A = 17^, and obtain the polynomial with real coefficients, 

P{X) = + c^X^ - C2X'^ - ciX + . (4.10) 

At a = 6 = the four eigenvalues of Ma,b,'y are known, which then gives 

P\a=b=o{X) = X^- 32(72 + 1)^2 + 256(1 - 272 + 7') . 

In addition, using the explicit formulas for the plane waves we compute the roots of P when 
a = 7 = 0: 

X^^'^^ = -4 ± 2v^6 + 562 ^ 0(^3)^ ^(3) ^ ^{A) ^ 4 _ 7^2 ^ 0(^3) _ 
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Similarly, when 6 = 7 = 0, we find 

= 4 ± 2V2a - 5a' + 0{a'), Xf) = Xj^^ = -4 + 7a' + 0{a') . 
Combining these formulas with the parity properties mentioned above, we conclude that 

C3 = 4(62 - a') + 0{a^ + 6^ + 7^), cs = 32 - 88(62 ^ ^2) ^ 33^2 ^ ^(^4 ^ ^4 ^ ^4) ^ 

ci = 0{a^ + 6^ + 7^), Co = 256 - 1664(62 + a') - 512j' + 0{a^ + 6^ + 7^) . 
A direct calculation now gives 

p(o) = 256 + 0(02 + 52 + ^2^ > , 

P(xf^) = -512a2 - 102472 + O((a2 + 72)(a2 + 62 + 72)) <0 , 

P(XW) = -51262 - 102472+ 0((62 + ^2)(^2^ ^2^ ^2)) ^0 ^ 

for a, 6, and 7 sufficiently small. This shows that the polynomial P has four real roots, so that the 
four eigenvalues of Aa^b,'^ are purely imaginary. This concludes the proof. □ 

4.5 Spectrum for 7 close to ^ 

In this case, the arguments are similar to the ones for (J2{A.a,b,'y) in Section 14.41 and we shall 
therefore only state the results and omit the proofs. First, we have the following result on the 
eigenvalues of A^^^ ^ . 

Lemma 4.9 There exist positive constants 71, c, and 5, such that the following hold, for any 
7 G [71,^], and {a,b) with ||(a, 6)|| < 6: 

(l) ia;+J^GB(-i;l),i^--lGi?(i;l); 

(^^) Kir i<^7 ' ''^tCf ^ ^(0; l),neZ\ {-1, 0}; 
(Hi) liw^;"^ - ^(^a,bj >c,n,peZ\ {0}, p ^ n, p -n; 
(^^) |i^<^C7 ~ 'Kai\ >c,n,peZ\ {-1}, p / n, p / -2 - n; 
M \^^aA^ - i<A7l > n e Z \ {0}, peZ\ {-1}. 
Next, we proceed as in the proof of Proposition 14.61 and obtain: 

Proposition 4.10 There exist positive constants 71, c, and 6, such that for any 7 € [71, and 
(a, 6) with II (a, 6)11 < 6, the following holds: 

(i) The spectrum of Aafirt satisfies 

a{AaM) C U B{\oJ-^l^^;c) U U B {iu;+l^^; c) , 

neZ neZ 

in which the balls B(\uj^'^^;c) are mutually disjoints except for pairs (\uj^'^^,\uj^'f^^) , n E 
Z\{0}, or(ia;-';^,iu;;-2-"), nGZ\{-l}. 
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(ii) Inside each ball B{iuj^'^^;c) the operator Aafi,^ has at most two eigenvalues, which are purely 
imaginary. 



5 The focusing NLS equation 

We consider in this section the focusing NLS equation 

\Ut{x,t) + U^^{x,t) + \U{x,t)\^U{x,t) = 0, (5.1) 

in which x S M, t G M, and f/(x, t) E C. This equation also possesses a family of small periodic 
waves of the form Ua,b{x) = e~'*e'^"'*'^Pa^b(A;a^;,x), but now 

4,6 = \{a'-h') + 0{a^ + h'') , 
ka,b = l + ^(a2 + 62) + o(a4 + 54) 
PaM = ae-'y + b^y + 0{\ab\{\a\ + \b\)) . 

Both the equation and the periodic waves have the same symmetry properties as in the defocusing 
case, so that we can investigate the stability of this family of periodic waves in an analogous way. 

As in Section 131 we define Pa^b and Qa,b{z) by (|3.1() . and consider solutions of 1)5. If) of the form 
(|3.3|) . The wave profile Qa,b{z) is then an equilibrium of the evolution equation 

iQt + ^Wa,bka,bQz + K,bQzz + (1 - plb)Q + \Q?Q = ■ (5.2) 

For the orbital stability, we use the same functional set-up, the same conserved quantities N{Q) 
and M{Q), and the energy 

r2-K 



£{Q) = ^ \2kl,\Q,{z)\'' - \\Q{zt) dz , 



in which only the sign of the last term has been changed. Following the arguments in Section |31 
one can show that the result in Theorem ^ holds in this case, as well. We only mention that the 
Hessian matrix of the function da b has now the expression: 



def / du)'^ divdc 
^ dc duj dcP 



Ul _\)(l + 0(a^+&^)) 



(a;,c)=(0,0) 

SO that it has again one positive and one negative eigenvalue. 

The analysis is also the same for the spectral stability, when we study the spectrum of the linear 
operator 

Aa,bQ = ^-^kl^bQzz - ^Pa,bka,bQz + i(l - pXb)Q + 2i|Qa,fe|^Q + "^Q^bQ ■ 

However, in this case the result is different: the periodic waves are spectrally unstable. While we 
do not attempt a complete description of the spectrum, we focus here on the existence of unstable 
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eigenvalues. It turns out that unstable eigenvalues arise through the unfolding of the quadruple 
zero eigenvalue of the unperturbed operator at a = 6 = 7 = 0. These are the eigenvalues of the 
4x4 matrix Ma,b,'y in Proposition 14.81 which is obtained here in the same way. This matrix has 
the same structure, 

_ /4i7(D2 + 0(a2 + 62)) M2ia,b) - A-f^l2 + Oia^ + b^)) 

but now 

The eigenvalues of Ma,b,'y are of the form A = i7X, where X is a root of a polynomial of the form 
(IHUj) . When a = 6 = we obtain 

P\a=b=oi^) = 32(7' + 1)^' + 256(1 - 272 + 7^) , 

and using the plane waves we find the roots of P when a = 7 = 0: 

'^^ = -4 ± 12^26 - 56^ + 0{h^), xf ^ = xj^^ = 4 + 76^ + 0{h^) , 

and when 6 = 7 = 0: 

= 4 ± i2^a + 50^ + O(a^), xf ^ = X^^^ = -4 - Ta^ + O(a^) . 

Then we find the expansions for the coefficients 

C3 = -4(62 _ ^2) ^ 0(^4 ^ 54 ^ ^4)^ C2 = 32 + 88(62 ^ ^2) ^ 32^2 ^ 0(^4 ^ ^4 ^ ^4)^ 

Cl = 0(a4 + 6^ + 7^), Co = 256 + 1664(62 + a^) - 51272 + 0{a* + b^ + 7^), 
which give 

P(0) = 256 + 0(02 + 62 + 72) > , 
Pixl"^) = 512a2 - 102472 + 0((a2 + 72)(a2 + 62 + 72)) , 

P{Xi^y) = 5 1 262 _ ;^Q24^2 ^ 0((^2 ^ ^2)(^2 ^ ^2 ^ ^2)) ^ 

This suggests that the polynomial P has complex roots provided 7 is small compared to a and 
6. In order to prove this, we consider the polynomial P when 7 = and show that it has at 
least two complex roots. We assume that 6 > a > 0, without loss of generality. Since P{X) = 
{X — 4)2(X + 4)2 + C)(a2 + 62), this polynomial is positive outside two C'(6-'^/2)_neighborhoods of 
4 and —4, when a and b are sufficiently small. Inside each of these neighborhoods, P has at most 
two real roots. A direct computation gives 

P(-4+y) = 51262 + (5l262+896a2)y+(64-4O62-l36a2)y2_(i6_4a2+462)y3+y4_^0^^4_^^4^) ^ 
from which we conclude that P is positive inside any 

0(6i/2) 

-neighborhood of —4, for ||(a, 6)|| 
sufficiently small. Summarizing, P has at most two real roots, and we conclude that the operator 
Aa,b,'y has at least one pair of eigenvalues off the imaginary axis, for 7 sufficiently small. In view of 
the symmetry with respect to the imaginary axis of the spectrum of Aa,b,'y, one of these eigenvalues 
has positive real part. This proves that the small periodic waves are spectrally unstable in this 
case. 
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Appendix: Spectrum of Ha^b 



In this Appendix we discuss the spectrum of the hnear self-adjoint operator Ha^b defined in 

As in Section |1J we decompose the elements of our function space into real and imaginary parts, 

and work with the matrix operator 

-4:kl bdzz + (Pa fe - 1) + 3i?a b + ^ab ^Pa,bka,bdz + 2Ra,bIa,b 



, _^^2 , 1\ I qd2 ,7-2 



-'^Pa,bka,bdz + 2Ra,bIa,b b^zz + {pI b ' + K b + ^/^ 



where Qa,b = Ra,b + ^Ia,b- We prove the following result: 

Proposition A.l There exists a positive constant sq such that for all {a,b) with ||(a,6)|| < eq, the 
spectrum of the matrix operator in the Hilbert space of 2tt -periodic functions Lp|,j.([0, 27r], C^) 
verifies 

a{Ha,b) = {0, a1% aI'J} U ai{Ha,b), <n{Ha,b) C [6,+oo) , 
where is a double eigenvalue and A^^, j = 2, 3, are simple real eigenvalues with 

>^Ofi = (^^d a2 < Ai% for all {a,b) + 0. 

Proof. Notice first that the parity properties with respect to (a, of the quantities kaf,-, Pa,b, 
and Qa^b imply that a{Hafi) = a{H„a,b) = (^{Ha,-b), and that Ha^b commutes with the symmetry 
S introduced in (|4.4j) . 

When a = b = the operator i/a,6 reduces to the operator Hq in the proof of Lemma 13.71 
with spectrum a{Ho) = {4n(?7- ±1), n G Z}, for which is a quadruple eigenvalue and the other 
eigenvalues are all positive and greater or equal to 8. Then, a standard perturbation argument 
shows that the spectrum of Ha^b decomposes as 

a{Ha,b) = {Xfl Ai^), X'fl} U ai{H,^b), where ai{Ha,b) C [6,+oo) , 

for (a, b) sufficiently small. The four eigenvalues A^"'! are the continuation for small (a, b) of the 
quadruple zero eigenvalue of Hq. 

In order to locate these four eigenvalues, we proceed as in the proof of Proposition 14.81 we 
construct an appropriate basis for the associated four-dimensional eigenspace, 

compute the 4 x 4-matrix Ma,b representing the action of Ha^b on this basis, and finally show that 
the eigenvalues of this matrix have the desired property. When a = 6 = we choose again the 
basis 

(0) _ / Sinz\ (1) _ /0\ (2) _ / COSZ \ (3) _ / 1 

For (a, b) ^ 0, the fact that Ha^b commutes with the symmetry S allows us to choose the vectors 
in the basis such that 

<;f(0) _ _.(0) ^^(1) _ Jl) ot(2) _ .(2) ^^(3) _ ^(3) 



and to conclude that the matrix A4a,b is of the form 

Mah = 



A2(a,6) 
B2(a,5) 
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where A2(a, b) and B2(a, b) are 2 x 2-matrices with coefficients of order 0{a'^ + 6^). 

Next, the two vectors dzQa,b and iQa.fe in the kernel of Aa,b also belong to the kernel of i^a,6, so 

(7) (7) (7) 
that we can take = q, for j = 0, 1, where q are as in from the proof of Proposition 14.81 

Then A2(a,6) = 0, so that zero is a double eigenvalue: A^'^^ = A^^^ = 0. The remaining vectors 

^ah j ~ ^' ^'^'-^ matrix 62(0, b) are computed from the expansions of ka^b, Pa,b, and Qa,b- We 

find ^2 = dfi + + b^) for J = 2, 3, and 

B2(a,6) = ^(3) (2), ,^ ,(3) ,(;)J+^(I«^K« +^)) 



\ \-"a,fe?o,0' ?0,0/ \-"a,fe?0,0' ?0,0/ 

2a^ Aab\ ,9xx 

4a* ^.^j+Oda^K.^ + i,^)). 



Since the spectrum of Ha^b is the same for all couples (±a, ±6), the determinant of 62(0, b) is even 
in both a and 6, which together with the above formula gives 

det (B2(a, b)) = -Ua^b'^ + 0(0^62(02 + 52^^ ^ q ^ 
This shows that A^^^ < < A^^'^^, for sufficiently small ||(a,6)||, which concludes the proof. □ 

Remark A. 2 We obtain the same result in the focusing case considered in Section\^ when the 
operator Ha b is given by 



Ha,b = £a,biQa,b) = -^K^b^zz - ^iPa,bKbdz - 0--Pa,b) " \Qa,b\ " 2Qa,fe ® Qa 

The only difference in the proof is the expression of the matrix 62(0, b) which is now 

B2ia,b)= -_ll^yO{\ab\{a' + b^)), 

but has the same determinant det (62(0, b)) = —I2a'^b'^ + 0{a%'^{a^ + 6^)) < 0. 
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